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To include the explicit valence gluon degrees of freedom into spectroscopy 
of the lowest states of heavy quarkonia, we consider, within an adiabatic 
model, the properties of the lowest hybrid QQg-mesons and estimate the 
effects of their mixing with low-lying vector cc-charmonia. The perspectives 
of compatibility of the resulting picture with data are discussed. 

The spectroscopy of systems constructed of a heavy quark Q and a heavy 
antiquark Q is known to be well described by the phenomenological potential 
models [|]J, where free parameters in the potentials are determined by fitting 
the calculated observables to the data. The gluon degrees of freedom are as- 
sumed to be integrated out in these models. The validity of this procedure is 
however questionable due to the presence of slow, large-scale nonperturbative 
vacuum fluctuations of the gluon field, as it follows from the lattice QCD ap- 
proaches. 

On the other hand, recent progress in understanding the production and de- 
cay processes of heavy quarkonia[|| is related to the idea of the valence gluon 
admixture in heavy quarkonium state vectors and the presence of the colour- 
octet QQ - subsystem inside the colour-singlet bound state 



Therefore, those quarkonium states which are outside the potential regime 
should not be used to fix parameters in the potential models, while for the 
states containing the gluon admixture the application of the potential models 
alone should leave a sufficient " room" for improvement assigned to subsequent 
inclusion of the gluon degrees of freedom. 

The key idea of our approach is nonperturbative mechanism of higher Fock - 
components ( i.e. the state vectors with the constituent gluons, Eq.(l) ) gen- 
eration in heavy quarkonia via the mixing with low- lying hybrid states. In the 
hybrid states, gluons are confined in a bound state by the confining interac- 
tion with the colour-octet quark core, which is assumed to be represented by 
an effective potential acting in three-body systems. Turning to this particular 
picture of the heavy hybrid composition, we notice existence of slow (QQ) 
and fast (g) sub-systems. Therefore, it is natural to proceed|| in the spirit 
of the Born-Oppenheimer, or adiabatic approximation, i.e. to solve first a 
relativistic wave equation for the gluon moving in a (presumably, confining) 
field of fixed center, and then to make use of the found gluon energy e g as a 
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part of the potential entering into the Schrodinger equation for the slow QQ 
- sub-system 



We note that the " Coulomb" potential is repulsive here because the QQ -pair 
is in the colour octet state. Our further main assumptions are as follows.We 
take the lowest magnetic (Ml)- and electric (El)- modes for the spin-orbital 
wave function of gluons 

^ = ^,m(^)b=l (3) 
= [^^J-l,m(0, <p) + |p4' + l,4», <f)]j=l (4) 

where Yjj, m (9, tp) are the vector spherical harmonics. This means that for the 
J PC = 1 - hybrid mesons we are going to consider, the orbital momentum I 
of a QQ-system should be I = 0(1) for the M1(E1) - gluon modes. The radial 
part of the gluon wave function is assumed to obey the Klein-Gordon (KG) 
equation with an external field including the (strong) Coulomb potential and 
the squared form of the linear confinement potential, properly scaled against 
analogous potentials for colourless QQ -states ( the scaling factor being the 
ratio of the corresponding Casimir operators equal to 9/4). For the assumed 
interaction between the two colour-octet, point-like particles, the KG-equation 
would be of the form 

(f g + V a ° 2 {r g ) + 2e g V v °{r g ) - V v ° 2 (r g ) - e g 2 )W 9 ) = (5) 

As far as the colour "charge" of the QQ - sub-system is spatially distributed, 
we define the "form factor-modified" potentials through the folding integral 

V s{v) (f g ) = J V s(v) (\ f g - r* |)p(rVr' (6) 

where the density function p(r) is related to an (unknown) wave function of 
heavy quarks. 

Finally, our calculation scheme acquires the variation form. We take, as a 
trial wave function of quarks, simple expressions of the exponential form ( 
with the pre-exponential centrifugal or nodal factors) which contain one vari- 
able parameter. This parameter propagates to the gluon energy, and then 
it appears again in the equation of motion for quarks. The last step is the 
minimization of the Schrodinger energy functional leading to the numerical 
value of this variable parameter and all energies,hence, to the hybrid meson 
mass. Summing up, with the trial wave functions of exponential form and on 
the basis of the adiabatic approximation, the masses of the vector hybrid 
states were estimated to be 4.02 (or 4.21) GeV for charmed quarks with mass 
m c =\A GeV and a valence gluon of the Ml(or El)-type,while for the bottom 
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quarks with mass mb =4.8 GeV the corresponding masses are 10.65 (or 10.75) 
GeV. The mean values of r g and Rqq characterizing spatial extension of hy- 
brid wave functions are equal to .45 fm and .4 fm for the charmed states, and 
.47 fm and .3 fm for the b-flavored hybrids. The approximate independence 
of the characteristics of light particle (i.e. the gluon) on masses of quarks 
is familiar manifestation of the heavy quark symmetry. The values of the 
"Coulomb" constant and slope of the linear potential have been taken equal 
to K=^a s =A9and&= .16 GeV 2 . The obtained values are rather close to esti- 
mates from different models §. In particular, they are very close to the string 
model elaborated in || where the gluon-quark interaction has acquired the 
form following from our formula (6) if 

p(r) = l/2(5{f-f Q ) + 5(f- r Q )) 

It is quite natural to expect that proper estimation of the hybrid meson(s) 
mixing with nearby quarkonia will be important to understand some pecu- 
liarities of the charmonium spectra and decays slightly over 4Gey||. 
As a first step, we confine ourselves to consideration of the four-level mixing in 
the charmonium spectrum choosing the ground state J/^(1S), ^(2S), ^(SS 1 ) 
and the presumed hybrid state H c with the calculated mass around 4GeV as 
mixing states. The nondiagonal elements m n H in the 4x4 - mass-matrix 

m nH = (H c ;QQg\TC int \nS;QQ) (7) 

i 

+ ^Y.\W(^- AM1 (r,e,m(r-m 
are calculated with the explicit radial wave functions 

R 3S (r) = N 3 ■ (1 - a 3 ( 7 r) m + b 3 ( 7 r) 2m ) exp (l/2( 7 r) m ), m = | (8) 

i?is = W3" -"r;a = & = 0), (9) 
#2S = W3" ^"2"; 6 = 0), (10) 

parametrized to reproduce approximately the spatial dimensions ( i.e.< r 2 >),the 
location of the radial function nodes and the values of the wave functions of 
the IS- 3S -quarkonia states at "zero" distance, which correspond to the 
QCD-motivited potentials (e.g. J?]] and references therein). 
As representative sets of [ 7n ; a n ; b n ] for the (nS)-states of charmonia we take 
[ 7l = .883; a x = h = 0], [72 = -715; a 2 = .57; b 2 = 0] [73 = .628; a 3 = 
1.101;63 = .197], where all 7's are in units of GeV. The lowest vector hy- 
brid state h c (gMiQQ) with the Ml-type gluon mode, i.e., with l g = 1 and 
Lqq = also called "the gluon-excited state", should presumably be rather 
narrow due to the dynamical selection rule discussed in a number of earlier 
papers^, |9|, which prevents the decays of this state into the ground-state 
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charmed mesons. This selection rule is not acting for hybrids with the El- 
type gluon mode (l g = 0, Lqq = 1), or "the quark-excited state" , which should 
therefore have very large width[|l0|. Hence, in what follows, we consider the 
mixing of the low-lying vector charmonia with the gluon-excited, Ml-type 
vector hybrid state. The radial wave function of this hybrid meson is taken 
in the factorized form in accord with the adopted adiabatic approximation 

RgQQ(r g ,R Q ) = N g N Q r g e W (l/2(a g r) m + l/2(f3 Q R Q ) m ),m = ~ (11) 

with the numerical values a g = 1.235 GeV, (3q = .973 GeV. The nondiag- 
onal elements of the symmetric 4x4- mass-matrix have been calculated as 
matrix elements of the interaction hamiltonian (7) over the reS'-charmonia 
states (n = 1,2,3), and the ("fourth") hybrid state and their values are: 
m n 4 = .25, .074, .044 GeV for n = 1,2, 3, respectively, with all other nondiag- 
onal elements equal to zero. 

To obtain " physical" eigenvalues of the diagonalized matrix near to the known 
masses of the ¥-family, we take the "bare" masses, which stand along the 
main diagonal having values Tn nn = 3.153; 3.695; 4.05 GeV for n = 1,2,3 
and m^{h c ) = 4.07GeV. The diagonalization procedure leads to the "phys- 
ical" masses: m(J/¥) = 3.089[3.097], m(¥(25)) = 3.685 [3.686], m(¥ (35)) = 
4.03(4.04], m(H c ) = 4.17(4.16] where masses of the known ¥-mesons in GeV 
are indicated in parentheses. The corresponding eigenfunctions reveal the 
following quark-gluon configuration mixing 

j/y = .968^(15) + .03020(25) + .0113^(35) - .2A7ip(h c ), (12) 

¥(25) = -.0628^(15) + .989^(25) + .0161^(35) - .134^(/i c ), (13) 

¥(35) = -.0908^(15) - .0697-0(25) + .940^(35) - .321ip(h c ), (14) 

H c = .223^(15) + .141^(25) + .332^(35) + .906ip(h c ). (15) 

Assuming the dynamical dominance of the quarkonia-components in the men- 
tioned states while computing the leptonic decay widths of the corresponding 
vector mesons, we can compare the model and phenomenological ratios of the 
meson wave functions "at the zero relative Q — Q-distance" : 

#W(0) = ^ls(O) : i?ls(0) : R 2 Hc (°) = 1 : - 65 [- 64 ] : -34[.27±.05] : .34[.29± .09], 

(16) 

where the values in parentheses are calculated using the proportionality be- 
tween Ry(Q) and m 2 (V)T(V — > The puzzling equality of the leptonic 
widths of the ¥(4.04) and ¥(4.16) states is explained in our approach by a 
coherent sum of the admixture (separately, looking small) amplitudes of the 
(15) — (35) quarkonia states in the dominantly hybrid ¥ (4. 16)-resonance. Whether 
our interpretation of this phenomenon can be experimentally distinguished 
from the Ono-Close-Page scenario^, |ll| (the approximately equal partition 
of the (35)-state and the H c -state between the ¥(4.04) and ¥(4.16) char- 
monium state vectors) remains to be considered. To deal with many other 
interesting implications of the valence gluon admixtures in the low-lying char- 
monia, the involvement problem of the broad (gEiQQ)-type hybrid state(s) 
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needs to be clarified. 

This work was supported in part by the Russian Foundation for Basic Re- 
search grants 96-02-18137 and 96-15-96423. 
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